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Abstract 

We study long-distance effects in rare exclusive semileptonic de- 
cays B — » K\£ + £~ , K\ is the axial vector meson. The form factors, 
describing the meson transition amplitudes of the effective Hamilto- 
nian, are calculated using Ward identites which are then used to calcu- 
late branching ratio and forward-backward asymmetry in these decay 
modes. The zero of forward-backward asymmetry is of special interest 
and provide us the precission test of Standard model. 

1 Introduction 

The investigation of rare semileptonic decays of the B meson induced by 
the flavor-changing-neutral-current (FCNC) transitions b — > s provide po- 
tentially stringent tests of standard model (SM) in flavor physics. In SM 
these FCNC transitions are not allowed at tree level but are induced by the 
Glashow-Iliopoulos-Miani (GIM) amplitudes pQ at the loop level. Addition- 
ally these are also suppressed in SM due to their dependence on the weak 
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mixing angles of the quark-flavor rotation matrix — the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [2J. These two circumstances make the FCNC de- 
cays relatively rare and hence important for the presence of new physics, 
commonly known as physics beyond SM. 

The experimental observation of inclusive [3] and exclusive [I] decays, 
B — > Xg'j and B —>■ K*^ , has prompted a lot of theoretical interest on rare 
B meson decays. However, in case of exclusive decays any reliable extraction 
of the perturbative (short-distance) effects encoded in the Wilson coefficients 
of the effective Hamiltonian [3 El [7J [HJ [9] requires an accurate separation of 
the nonperturbative (long-distance contributions), which therefore should be 
known with high accuracy. The theoretical investigation of these contribu- 
tions encounters the problem of describing hadron structure, which provides 
the main uncertainty in the predictions of exclusive rare decays. In exclu- 
sive B — ► K, K* decays the long-distance effects in the meson transition 
amplitude of the effective Hamiltonian are encoded in the meson transition 
form factors. Many exclusive B -> K (K*) £ + £~ [TOJ QH E2] , B -> 7^" [13], 
B — > £ + £~ [H] processes based on b — > s (d) £ + £~ have been studied in liter- 
ature and many frameworks have been applied to the description of meson 
transition form factors: among them the worth mentioning are constituent 
quark models, QCD sum rules, lattice QCD, approaches based on heavy 
quark symmetry and analytical constraints. Many observables like Forward 
Backward (FB) asymmetry, single and double lepton polarization asymme- 
tries associated with the final state leptons, have been extensively studied 
for quite some time for quark level processes b — > s (d) £ + £~ . 

Recently, Belle[15] has announced the first measurement of B — > K± (1270)7 

B(B + -> K+i) = (4.28 ± 0.94 ± 0.43) x 10~ 5 . (1) 

after which these radiative decays became topic of prime interest and their is 
lot of theoretical progress for which one can see the Refs. fTBlfTT] . In this paper 
we study the semileptonic B meson decay B — ► Ki£ + £~ using the framework 
of Gilani et al. [IB] with K x is an axial vector meson. The axial vector mesons 
is distinguished by vector by an extra 7 5 in the gamma structure of decay 
amplitude (DA) and some non perturbative parameters. But the presence 
of extra 7 5 does not alter the calculation except the switching of vector 
to axial vector form factors and vice a versa. As mentioned earlier, the 
theoretical understanding of exclusive decays is complicated mainly due to 
non-perturbative form factors entered in the long distance non-perturbative 
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contributions. The aim of this work is to relate the various form factors in 
model independent way through Ward identities. This enables us to make a 
clear separation between non-pole and pole type contributions, the q 2 — > 
behavior of the former is known in terms of a universal function £ ± (0) 
tj . (0) introduced in the large energy effective theory (LEET) of heavy (B) 
to light (Ki) form factors[T7]. The residue of the pole is then determined in 
a self consistent way in terms of g + (0) or £ ± (0) which will give information 
about the couplings of B* (1~) and B* A (1 + ) with BK\ channel. The from 
factors are then determined in terms of the known parameters like g + (0) 
and the masses of the particles involved which are then used to calculate the 
branching ratio and forward-backward asymmetry for these decays. 

This paper is organized as follows: In section II we introduce the effective 
Hamiltonian formalism of semileptonic B meson decays and will write down 
the matrix elements for B — > K\£ + £~ decays. Section III discusses the Ward 
identities and develop the relationship between form factors which results 
in the reduction of number of unknown quantities. The form factors thus 
obtained are used for the calculation of decay width and forward-backward 
asymmetry. Finally, in the last section we summarize our conclusions. 



2 Effective Hamiltonian 

At quark level the decay B — ► Ki£ + £~ is similar to one studied in, for 
example, reference [10]. The basic transition b — > s£ + £~ is described by 
the effective Hamiltonian given below 

c 10 

H eff = -4-| V th V* J2 CiMOiO*) (2) 

* i=l 

where O^s are four local quark operators and Cj are Wilson coefficients 
calculated in Naive dimensional regularization (NDR) scheme |19j . 

One can write the above Hamiltonian in the following free quark decay 
amplitude 

c e 9 ff [s lfl Lb] [hn] ) 

+C 10 [s llx Lb] [i^H] \ 
2m b C e 7 ff [sia^Rb] [t^t] J 

(3) 



M{b 



s£ + t 



Gpa v V* 
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with L/R= s = q 2 which is just the momentum transfer from heavy 

to light meson. The amplitude given in Eq. ([3]) is a free quark decay amplitude 
which contains certain long distance effect from the matrix element of local 
quark operators, (l + l~s |Oj| b) , 1 < i < 6, which usually reabsorbed into the 
redefinition of short distance Wilson coefficients. Specifically, for the exclu- 
sive decays, the effective coefficients of the operator Og = -^^(sj^Lb)^^) 
can be written as 

Cl fS = C 9 + Y(s) (4) 
where the perturbatively calculated result of Y(s) is [T9l 120] 

g (m c ,s) (3Ci + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) 
^ pert (s) = -\g (1,§) (4C 3 + 4C 4 + 3C 5 + C 6 ) (5) 

-f <? (0,s) (C 3 + 3C 4 ) + | (3C 3 + C 4 + 3C 5 + C 6 ) . 

For the values of the Wilson coefficients and the explicit expressions of g's 
appearing in Eq. ([5]) we will refer to [HI [20] . The hat denotes normalization 
in terms of the B meson mass[10j. 



3 Matix Elements and Ward Identities 

Exclusive decays B — > Ki£ + £ involve the hadronic matrix elements of the 
quark operators in Eq. ([3]) between B and K\. These can be parameterized 
in terms of form factors which are the scalar functions of the four momentum 
square (q 2 = {j>b — PrJ 2 )- For the process we are considering, there are 
seven form factors like the transition of psudoscalar to vector meson. The 
non-vanishing matrix elements are 

(K 1 (k,e)\V,\B(p)) = ie;(M B + M Kl )V 1 (s) 

-{P + k),{e*-q) V - [s) 



M B + M Kl 



-%{e-q) 2 -^ms)-V {s)] (6) 
(K 1 (k,e)\A,\B( P )) = J ie ^ e-p^A(s) (7) 

1V1 B + M Kl 

with = sj^b and = s^^^b are the vector and axial vector currents 
respectively and e* is the polarization vector for the final state axial vector 
meson. In Eq.dH]) 
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V*(s) 



M B + M K M B -M Kl jrr 



2M, 



K, 



2M, 



(8) 



with 



V 3 (0) = V Q (0). 

In addition to the above form factors there are also some penguin form factors 
which are: 

(K l (k,e)\sia, v q»b\B{p)) = [(M% - M|J e M - (e ■ q){p + fc)J F 2 (s) 



He* ■ q) 



M B -M 



-(p + fc), 



(10) 



(K 1 (k,e)\sia^j 5 b\B(p)) = 
with 

F 1 (0) = 2F 2 (0). 

The various form factors appearing in Eqs. (|6l)-(|T0l) can be related by Ward 



-{m b + m s )(K 1 (k,e)\sj ti b\B(p)) (11) 
m b - m s ) (K x (h,e) \si^ 5 b\ B{p)) 



identities as follows [T8| 

(Ki(k, e) \sia flu q u b\ B{p)) = 
(Kx(k,s) \M<r /a ,q u 'y 5 b\ B(p)) = 

+(p + k) fl {K 1 (k,e)\sj 5 b\B(p)). (12) 

Now we make the heavy quark approximation and compare coefficients of £* 
and g M from both sides. In the heavy quark approximation we need not to 
compare the coefficients (p + k)^. Using Eqs. (|6l)-(|T0l) in Eqs. (TTTT) and ( |T2l) . 
we get the following relationship between form factors 

[m b - m s ) 



F 2 (s) 
F«(s) 



M B + M Kl 
(m b + m 



2A(s) 
Vi(s) 



M B - M Kl 

(m b + m s ) [V 3 (s)-V (s)\. 



2M Kl 



(13) 
(14) 
(15) 



These are model independent results derived by using Ward identities. The 



universal normalization of the above form factors at q 2 
by defining [T8] 

(K 1 (k 1 e)\sta a/3 ^b\B(p)) = 



-16 



a/3pa * 



is obtained 



< [(P + k) p9+ + q a g_] - (q ■ e*)e^{p + k) p q a h 



-i [q a e^s;(p + k) a q T 



a: 



(3] h. 



(16) 



Using Dirac identity 



(17) 



in Eq. ([16]) one can write 



(ifi(fe,e) -B(p)) = ^[(Ml-A/IJ (/+ + */_] -(g-0[(P + fcW + ?^-] 

+(g • e*) [a(p + - (Mi - M^JgJ /i. (18) 

Comparing coefficients of g M , e* and e Ml/a/ g from Eqs.flH]), f TTUT) and Eqs. ( |T6i) 
and (|T8|) . we get 



Fi(s) 



2 [g+(s 

9^ 



s 



Ml 



Ml 



-g. - {M 2 B - M 2 Kl ) h. 



(19) 
(20) 

(21) 



The above results ensure that iq(0) = 2F 2 (0). In terms of g + , g>_ and h, the 
form factors become 



A(s) 
Vi(s) 
V2(s) 



M B + M Kl 
m b - m s 
M B -M. 



\g+(s) - s/ii] 



A'i 



m b + m s 
M B + M Kl 
m b + m s 



9+ + 



Ml 



Ml 



[g + (s) - sh] 



2M 



Ki 



Mr — M, 



-V (s). (22) 



By looking at the above expressions one can see that the normalization of 
above form factors A and V± at s — is determined by the single constant 
g + (0) where as that of V 2 is determined by g + (0) and V (s). 



3.1 Pole contributions 

The pole contribution for B to p has been studied in detail by Gilani et 
al. [18]. This remains the same for B to K\ transition and again only hi, 
<7_, h and Vq get pole contributions from B*(l~), B* A {1 + ) and B(0~) mesons 
where as g + , g- and Vq(s) gets their contribution from quark triangle graph. 
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These are given by 

hl\pole 
9— \pole 

h\pole 



B* 
T 



2 M|* 1 - s/M 2 B . M|, 1 - s/M%, 
9b* a bk x f 



R A 



M 2 B , 1 - s/M%* M 2 B , 1 - s/M 



1 Jb*. BK\ 



It 



R L 



2 M 2 * l-s/M B * 
V (s)\ po u 



M 2 1 



Qbbk\ f s/M B 

JB 



Rr] 



s/M 2 



M Kl JD 1- S /M 2 " u 1 - s/M! (23) 

where Ry, R A , Ra an d -^o are related to the coupling constants 9b*bk^-, 
9b\bk y -, f b* a BK\ and 9bbk x respectively. One can find the detail about it in 



Ref.|18j. Thus one can write 
A(s) > M * + M ^ 

Vi(s) 



m b - m s 

M B - M Kl 

m;, + m s 



V2{S) 



M B + M Kl 
m b + m s 




The behavior of g+(s), g~(s) and Vq(s) near s — > is known from LEET and 
their form is 1 1 8 1 



g+(s) 

V (s) 



e±(o) 



2\2 



(1 - s/Ml) 

1 



M 2 Kl 



At s -> 



Vb(0) 
9+{0) 



M B E Kl 



Ml -M 2 K 



-9-{ s ) 



M 



2 

e±(o). 



Ml 



+ 



(27) 
(28) 

(29) 

(30) 
(31) 
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The pole terms in the relations (1241) . (T2"5"l) and (1261) are expected to dominate 
near s = or M^* . On the other hand the relations obtained from Ward 
identities, are expected to hold for s much below the resonance region. The 
above behavior, near s = and that near the pole [18] suggest 

F(s) = m m 

y } (1 - s/M 2 ) (1 - s/M' 2 ) y } 

where M 2 is M\* or M§* and M' is the radial excitation of M. This pa- 
rameterization not only takes into account the corrections to the single pole 
dominance, as suggested by dispersion relation [22j [231 [21], but also of off- 
mass-shell-ness of couplings of B* or B* A with BKi channel. 

Since g+(s) and g~(s) have no pole at s = M^„, therefore we get 

V MJ,« J B * \ m b -m s 

This gives using the parametrization fl32l) 

Rv = -~9b*bkJF = -\g B *BKj B = - 1 _M2° t/M a ( 33 ) 



Similarly, 

R A = o f b*. BK%fx^ = —, 1*2 '1**12 ■ (3 4 ) 



,;»^^ f B* A _ <y+(0) 



2 JB* A BK ljT 1-M 2 */M' 2 ,' 



For the detailed derivation and discussion on these relations we will refer to 
[T8] . We cannot use the parametrization given in Eq. (1321) for Vi(s) since near 
s = 0, 7i(s) behaves as g+(s) [l - s/(M| - M 2 Ki )] [c.f. Eqs. (ESI) and ([271)]. 
This suggests the following 

VM = 7 ( 1 - ,,^^1 . (35) 



(l - S /*fi.) (l - s/A/g. ) V Ml - Mf, 

Until now we have expressed every thing in terms of g+(0) which is the only 
unknown in the calculation. After the first announcement of Belle [15] for the 
decay B — > Kij, the value of g+(0) has been extracted to be [HI [TF] 

9+(0) = S±(0) = 0-32 ±0.03. (36) 
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Using Jb = 180 MeV we have the prediction from Eq. 

g B *B Kl = 15.42 GeV- 1 (37) 
Similarly, the S and D wave couplings are predicted to be 

9b a bk x = S.17f B1BKl GeV 2 (38) 
The different values of F(0)'s are 

A/ . (M B + M K , \ 

A (0) = _? -)g+(0) (39) 



mi, — m 



s 



: m b - m k . , 

^i(O) = - 1 ■ * 9+(0) (40) 



rrib + 7n 



s 



*(„) = _ (Mb + MkA (q) _ _^!Ll_v (0) (41) 
\ m b + m s J M B - M Kl 

where <?+(0) is the same as given in Eq. ( 1361) . The calculation of numerical 
values of A(0), Vi(0) is very trivial but to go for V^O), we have to know the 
value of Vq (0). Although LEET does not give any relationship between £||(0) 
and £j_(0) but due to some numerical coincidence in the LCSR expressions 
for f,|(0) and ^(0)[2S] 

^ n (°)-M0) = P+(0) (42) 

so that, from Eq. (J2HD 

V (0) = 1.13^(0). (43) 

Thus the final expressions of form factors which we shall use for numerical 
work are 

Mo) 



M*) 



s/M%) (1 - s/M'i) 



Vi(s) = -7 ^ rx(l-TT^T^) (M) 



Ml - Ml 



V 2 (s) 
where 



(l - S /Mi») (l - a/M§. 

V 2 (Q) 2M Kl V Q (0) 



(l - s/M B1 ) (l - a/M^j M b ~ M k, (1 - s/M B ) (1 - s/M'i) 

A(0) = (0.52 ±0.05) 
Vi(0) = -(0.24 ±0.02) 

V 2 (0) = -(0.39 ±0.03) (45) 
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4 Decay Distribution and Forward- Backward 
Asymmetry 

In this section we define the decay rate distribution which we shall use for 
the phenomenological analysis . Following the notation from ref-JTU] we can 
write from Eq. Q 



M 



GpOL 



v tb v*m B [t; (hn) + r M 2 (zWOl 



(46) 



where 

Tl = A {§) e wa ^ p v%VK, ~ iB (§) e* + iC (s) (e* ■ p B ) p hfl + iD (s) (e* ■ p B ) g M 

(47) 

CT 2 = E(s)e» pa ^p B f Kl -iF(s)e; + iG(s)(e*-p B )^ 

(48) 

The definition of different momenta involved are defined in reference |10j. 
where the auxiliary functions are 

■Cf{s) + 2 ^CfF 1 {^) 



A(s) = 


2A(s) 


l + M K 


B(s) = 


(l + M Kl 


C(s) = 


1 



Cl f \s)Vi(s) + —C e 7 ff (l-M Kl 



C e 9 ff {s)V 2 (s) + 2m b G 



l-M 2 Kl 



eff 



F 3 (s) + 



1 — Ml 



Fo(s) 



E(s) 

F(s) 
G(s) 

H(s) 



Cl ff (s)(l + M Kl )Vi(s) - (1 - M Kl )V 2 (s) - 2M Kl V (s] 



-2m b C e 7 ff F 3 (s) 



2A{s) 



-C 
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l + M Kl 
1 + M Kl ) CMs) 



-C 10 V 2 {s) 



l + M Kl 

I \c 10 (s)(l + M K i)Vi(§) - (1 - M Kl )V 2 (s) - 2M Kl V Q (s)] . (49) 
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The differential decay rate for B — > K*jj, + /j,^ can be expressed in terms of 
these auxiliary functions in[10j and this remains the same for B — > Kifi + {i~ 
with the obvious replacements. Integration on s in the range 

(2m) 2 < s < (1 - m Kl f (50) 

with rhi = mi/nriB, and using Tb° = (1-530 ± 0.009) x 10~ 12 s, the branching 
ratio is 

B (B -f Kifi + fi-) = 0.9±g;ii x 10~ 7 

The above value of branching ratio is for the case if we do not include Y(s) in 
Eq. (j3J). The error in the value reflects the uncertainty from the form factors, 
and due to the variation of input parameters like CKM matrix elements, 
decay constant of B meson and masses as defined in Table I. 

Table I : Default value of input parameters used in the calculation 



m w 


80.41 GeV 


m z 


91.1867 GeV 


sin 2 6w 


0.2233 


m c 


1.4 GeV 


Wlb,pole 


4.8 ± 0.2 GeV 


m t 


173.8 ± 5.0 GeV 


a s (m z ) 


0.119 ±0.0058 


fs 


(200 ± 30) MeV 


\v; s v tb \ 


0.0385 



Now if we include the value of Y(s) the central value of branching ratio 
reduces to 

B (B -> Kxn+n-) = 0.72 x 10~ 7 

By including Y(s) the behavior of the differential decay rate as a function of s 
is shown in Fig. 1. The solid line denotes the theoretical prediction with input 
parameters taken at their central values, while the band between two dashed 
line shows the uncertainity from input parameters. In our numerical analysis 
we have considered only the final state leptons as being the muon. Our reason 
for choosing this is due to the extreme difficulty in detecting electron in the 
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final state and that the branching ratio B — > K^i^ becoming small with 
the SM for the r in the final state. 

The differential forward-backward asymmetry for B — ► K\[i + reads as 
follows [TO] 

gL4 fb G 2 F a 2 m% 2 



ds 2 10 7T 5 

where 



V* s V tb \ su (s) [Re (BE*) + Re (AF*)] (51) 




u(s 



- 2m 2 Ki (1 + s) (52) 

The variable u corresponds to 6, the angle between the momentum of the B 
meson and the positively charged lepton in the dilepton cm. system frame. 
The behavior of forward-backward asymmetry in B — > Kifi + fi~ decay as 
a function of s is shown in Fig. 2. Contrary to the branching ratio, the 
forward-backward asymmetry is less sensitive to the input parameters as is 
clear from Fig. 2. For the zero-point of forward-backward asymmetry in the 
standard model, we get s = (0.16 + 0.01) (s = (4.46 + 0.27) GeV" 2 ). 
Conclusions 

We have studied B — > Ki£ + £~ decay using Ward identities. The form 
factors have been calculated and found that their normalization is essentially 
determined by single constant g + (0) which has the value g+(Q) = 0.32 ±0.03 
obtained from [T61 ITT] . By considering the radial excitation of M (where 
M = Mb* or M B * A ), which are suggested by dispersion relation [Tg]. we have 
predicted the coupling of B* or B* A with BK\ channel as indicated in Eq. 
( ]3Tj) and the value is gg* bk\ = 15.42 GeV -1 . Also we have predicted the 
relationship between S and D wave couplings gB* A BK± = 3.17/B^BE'iGeV 2 
given in Eq. (1381) . We have summarized our form factors in Eq. (14*41) and 
their value at s = in Eq. (I45p . By using these form factors we have 
calculated the branching ratio for B — > K\y^ both by considering the non 
resonant and resonant value of the Wilson coefficient C$ (s) which will been 
seen in future experiments. The decay distribution is shown graphically in 
Fig. 1, where the differential decay rate is plotted as a function s. 

A detailed analysis of the forward-backward asymmetry is also presented 
here. We have plotted the forward-backward asymmetry as a function of s 



12 



in Fig. 2. It is clear from the graph that the the SM the central value of the 
zero of the FB asymmetry is at s = 0.16 (s = 4.46). This value of the zero 
of the forward-backward asymmetry will provide the precision test of SM in 
planned future experiments. 
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Figure Captions 

1): The differential decay rate as a function of s is plotted using the form 
factors calculated by using Ward Identities. The resonanant cc states are 
parameterized as in refs. [T9l [20] . Here the solid line denotes the theoretical 
predictions with the input parameters taken at their central values, while the 
dashed (dotted) line is for max. (min) value of input parameters. 
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2): The forward-backward (FB) asymmetry as a function of s is plotted 
using the form factors calculated by using Ward Identities. The resonanant 
cc states are parameterized as in refs.[19, 20J. The dashed (solid) line is for 
the central (max.) value of the input paramteres. 
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